Decoding mode-mixing in black-hole merger ringdown 
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Optimal extraction of information from gravitational-wave observations of binary black-liole coa- 
lescences requires detailed knowledge of the waveforms. Current approaches for representing wave- 
form information are based on spin-weighted spherical harmonic decomposition. Higher-order har- 
monic modes carrying a few percent of the total power output near merger can supply information 
critical to determining intrinsic and extrinsic parameters of the binary. One obstacle to constructing 
a full multi-mode template of merger waveforms is the apparently complicated behavior of some of 
these modes; instead of settling down to a simple quasinormal frequency with decaying amplitude, 
some \m\ ^ £ modes show periodic bumps characteristic of mode-mixing. We analyze the strongest 
of these modes - the anomalous (3, 2) harmonic mode - measured in a set of binary black-hole 
merger waveform simulations, and show that to leading order, they are due to a mismatch be- 
tween the spherical harmonic basis used for extraction in 3D numerical relativity simulations, and 
the spheroidal harmonics adapted to the perturbation theory of Kerr black holes. Other causes of 
mode-mixing arising from gauge ambiguities and physical properties of the quasinormal ringdown 
modes are also considered and found to be small for the waveforms studied here. 

PACS numbers: 


I. INTRODUCTION 

Since the first successful simulation of black-hole bina- 
ries (BHBs) through late inspiral, merger, and ringdown 
in 2005 [1-3] , theoretical interest has centered on the re- 
sulting gravitational waveforms. A crucial tool in wave- 
form studies has been the analysis of the radiation wave 
pattern in spherical harmonic components. This decom- 
position is useful both in the physical interpretation of 
the radiation, and in structuring the waveform informa- 
tion content for the development of approximate analytic 
or empirical encodings. 

The self-consistency of results for the dominant 
quadrupole waveforms across numerical codes was 
quickly established [4, 5], enabling rapid study of the 
basic characteristics of mergers [6-13] Researchers soon 
began to build analytic template models compatible 
with these numerical results as well as with the post- 
Newtonian (PN) at earlier times, to provide relatively 
quick waveforms for specified BHB source masses and 
spins [14-16]. While expected to be sufficient for detec- 
tion of BHB mergers, quadrupole-only templates will not 
lock down most of the intrinsic (masses, spin magnitudes 
& directions) and extrinsic (sky position, phase) BHB 
system parameters. To gain an understanding of these 
parameters requires a richer template bank, one that in- 
cludes all of the relevant angular modes of the signal 
[17-20], 

Working with a spherical harmonic basis of spin- weight 
s = —2 [21, 22], several studies [23-27] have found that 
after the dominant quadrupole (£ = 2, m = ±2) modes, 
the next most important modes tend to be the higher 
m = ±£ modes: (3, ±3), (4, ±4), etc., though odd-m 
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FIG. 1: Gravitational-wave luminosity from the merger and 
ringdown of a nonspinning black- hole binary of mass ratio 4:1, 
decomposed by harmonic mode. 


modes are sometimes suppressed by symmetry. We have 
also seen, however, that certain m < \l\ modes can be 
important. Prominent amongst these are the (2, ±1) and 
(3, ±2) modes. Figure 1 shows the radiative power for 
the most important modes in the case of the merger of a 
4:1 nonspinning BHB. Here we see that the (2, 1) mode 
has actually overtaken the (5, 5) mode in importance by 
merger time. 

A key feature of BHB mergers exposed through the 
spherical harmonic decomposition waveform studies is 
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the rather clean separation of the sometimes compli- 
cated mix of signal frequencies, achieved by angular 
mode decomposition. Even when typical observers would 
measure complicated wave shapes combining several fre- 
quency harmonics, these harmonics largely reduce to 
slowly evolving sinusoids in each spherical harmonic com- 
ponent mode. To a very good approximation, this struc- 
ture holds consistently through the inspiral, merger, and 
ringdown [23, 28-30]. This pattern of frequency separa- 
tion is extremely convenient in allowing relatively simple 
encodings of the waveform information in analytic mod- 
els. 

Partly because of these properties, angular-mode de- 
composition has become a standard approach to compar- 
ing waveform simulations with each other, with analytic 
post-Newtonian calculations, and with developing empir- 
ical waveform template models. These uses of the decom- 
position technique have elevated its significance from its 
beginning as an interpretive convenience to its current 
status as an essential component of how we quantita- 
tively understand gravitational-wave signals. Thus we 
must be aware of the possibility that artifacts of arbi- 
trary choices in the details of the decomposition proce- 
dure may interfere with our quantitative understanding 
of the waveforms themselves. 

Such concerns are particularly notable when we see un- 
usual features in the decomposed waveforms seeming to 
violate the a posteriori expectation of clean separation 
of frequencies. Several authors [7, 25, 26] [BJK: who 
else?] have noted that the (3, 2) mode in partic- 
ular typically seems to break from this simple pattern, 
showing unusual post-merger features that require inves- 
tigation and resolution before a useful model can be de- 
veloped. In some of the earliest merger simulations, Buo- 
nanno et al. [7] already noted the presence in the post- 
merger “ringdown” (3,2) mode of both (3,2) and (2,2) 
quasinormal-mode (QNM) frequencies. 

In this paper, we investigate these (3, 2)-mode anoma- 
lies, with a survey of 3D numerical simulations of the 
merger of various comparable-mass BHBs with non- 
precessing spins, exploring a range of possible “causes” . 
We find that the dominant part of the measured mode- 
mixing that underlies the anomalous effect can be at- 
tributed to our use of spherical harmonics rather than 
the spheroidal harmonics expected by Teukolsky pertur- 
bation theory. 

The remainder of this paper is laid out as follows: In 
Section II, we review the numerical evidence for mode- 
mixing in existing (3, 2) evolutions, and show how well 
it is captured by a simple two-mode phenomenologi- 
cal model for the ringdown waveform segment. In Sec- 
tion III, we discuss general models for why mode-mixing 
should be expected, including effects of coordinate dis- 
tortions in the radiation extraction spheres, and of ill- 
adapted harmonic basis functions in the radiation decom- 
position. In Section IV, we introduce our set of expanded 
numerical evolutions, arranged into “equivalence classes” 
of common end-state Kerr spins, which we analyze in Sec- 


tion V, fitting the measured contributions of two-mode 
models to our models. We conclude in Section VI with 
discussion on the application of these results to more gen- 
eral late-merger-ringdown models, such as the implicit 
rotating source model of Refs. [25, 26]. We present a de- 
tailed description of our selection of equivalence classes 
of binaries in B. 


II. BUMPS IN NUMERICAL (3,2) MODES 

The first gravitational waveforms extracted from nu- 
merical simulations were the dominant (2, ±2) modes, 
whose early-inspiral behavior was expected to match the 
quadrupole radiation predicted by quasi-Newtonian and 
post-Newtonian theory. Once these had been shown to 
be robust and universal across codes [4, 5], some groups 
turned their attention to the subdominant modes. When 
Baker et al. [25] looked at a set of mergers of nonspin- 
ning black- hole binaries with mass ratios in the range 1:1 
to 6:1, they noted that one of the leading subdominant 
modes, (3,2), showed an unusual bumpiness just after 
merger over a range of parameter space. This bumpi- 
ness manifested in both the frequency and amplitude, 
and appeared to persist with both increased resolution 
and extraction radius, thus constituting a robust pattern 
of excursions from the frequency separation dominating 
the £ = m modes. More recent work by Kelly et al. [26] 
shows the same anomaly in equal-mass binaries with non- 
precessing spins (i.e., the spins are aligned/anti-aligned 
with the orbital angular momentum). Meanwhile, Buo- 
nanno et al. [7] report that an accurate fit of the (3, 2) 
mode for the ringdown stage of effective-one-body (EOB) 
waveforms requires the addition of the fundamental (2,2) 
quasinormal frequency. 

Examples of these more complicated waveform fea- 
tures are shown in Fig. 2, where we plot waveform fre- 
quency (top panel) and amplitude (bottom panel) of the 
measured (3, 2) mode for the merger of a nonspinning 
4:1 binary, as well as for the mergers of several other 
BH configurations with the same final dimensionless spin 
(«f fa 0.475). We also mark the expected real QNM (2, 2) 
and (3, 2) frequencies, W22 and W32 for a Kerr black hole 
of this spin. From the time of peak amplitude (f = 0 
here) until the waveforms start to degrade around 60M 
later, the frequency seems to oscillate around one or other 
of these two QNM frequencies, rather than locking onto 
the higher W32, as for other modes. These oscillations ap- 
pear in the strain h and its time-derivatives; we choose to 
study strain-rate, h(t), waveforms, which we decompose 
into modes fi(£, m )(t) = |h(f,m)| exp(i(^^ jTn )), with instan- 
taneous frequencies 

We can model the more complicated ringdown wave- 
form features by expressing the (2, 2) mode as a pure 
QNM ringdown, and the measured (3, 2) mode as a lin- 


3 



FIG. 2: Post-merger frequency [top panel] and amplitude 
[bottom panel] of the numerically measured (3, 2) mode for 
a set of “4:l-equivalent” evolutions, resulting in a final black 
hole with dimensionless spin at « 0.475, matching that of 
a 4:1 nonspinning binary merger. The data sets have been 
shifted in time so that t = 0 corresponds to peak amplitude 
of the dominant (2, 2) mode. The two dashed (black) hor- 
izontal lines in the top panel mark the fundamental QNM 
frequencies U 22 (lower) and UJ 32 (higher) for a Kerr hole of 
the same final spin. 

ear combination of QNM ringdowns: 

h$°. 2 d ) el = A 22 e i{ ° 22t+5 ™\ (1) 

h^°. ff = A 32 e^ t+5 ^ + p 32 A 22 e l ^ t+5 ^ . ( 2 ) 

Here = u!( m + */ r ^m ' s the full complex QNM fre- 
quency, and /J 32 = poexp(*C) is a constant complex- 
valued parameter indicating the mixing of the (2, 2) QNM 
mode into the measured (3, 2) mode. The modeled (3, 2) 
mode frequency and amplitude are then: 

, ■..model/.*', 1 A R f o\ 

(3,2) W - W 32 H (3) 

e(f) [A R cos(A R t + 5) + Aj sin( A R t + 5)] 

m 

|hf 3 ? 2 d ) el (t)| (4) 


where F(t) = 1 + 2 s(t) cos( A R t + 5) + e(t) 2 , A R = uj 32 — 
w 2 2 , A/ = l/r 32 - l/r 2 2 , e(t) = pqA 22 / A 32 exp(A/t) = 
£ 0 exp(A/t), and 5 = 5 32 — 5 22 — (. 

For a given mass and spin, the QNM frequencies, W 22 
and W 32 , and the damping times, T 22 and T 32 , are values 
known from black hole perturbation theory. Typically, 
T 22 « T 32 , so that Aj is somewhat smaller than A R , 
allowing a beat-like effect to persist over several cycles. 
Fixing these leaves just two free parameters for the fre- 
quency: £0 = PqA 22 / A 32 , the initial ratio of contributing 
amplitudes, and 5 , the initial phase difference, as well as 
one more amplitude parameter, A 32 . 

Evidently, the characteristic shape of the modeled (3,2) 
mode frequency plots will depend on the relative magni- 
tude of the modal contributions: for £ 0 1 , the fre- 

quency will oscillate (approximately) sinusoidally about 
W 32 ; for £q 1 , the oscillation will be about ui 22 ] for 
intermediate values, the oscillatory shape will be more 
complex. In the left panel of Fig. 3, we demonstrate 
these shapes for a Kerr hole of spin a{ = 0.475, the same 
4:1 end-state spin as in Fig. 2. Similarly, the right panel 
shows the corresponding modal amplitude shape for the 
same end-state hole. Again, the most extreme bumps 
in amplitude occur when the (2, 2) and (3, 2) modes have 
comparable amplitude contributions (£0 ~ 1). These the- 
oretical curves should be compared with the numerically 
measured mixing in Fig. 2. 


III. POSSIBLE CAUSES OF MODE-MIXING 

The bumpy features seen in the measured (3, 2) mode 
are a clear exception to the general rule that each angu- 
lar mode encodes sinusoidal waves with just one slowly 
evolving frequency component, the phenomenon we re- 
fer to as frequency separation. In Fig. 3, we showed 
that a combination of the fundamental (3, 2) and (2, 2) 
quasinormal-mode frequencies produces similar features. 
More generally, there are indications that such mixing 
occurs among other modes, especially other higher-order 
to = 2 modes, which likewise seem prone to coupling to 
the dominant mode. Here we ask the basic question: is 
this mode-mixing a fundamental property of the radia- 
tion, or some kind of an artifact, and if so, what kind? 

We consider various hypotheses to explain this mode- 
mixing effect violating our empirical frequency separation 
rule. The first, which we label physical mixing , is simply 
that the frequency separation rule does not physically 
hold to sufficiently high precision; that is we are are per- 
haps seeing a nonlinear effect in the radiation-generation 
process underlying the (3, 2) mode. Under this assump- 
tion, no choice of fixed or slowly evolving angular basis 
could be expected to yield the kind of frequency sep- 
aration we see in other cases. Near the merger where 
nonlinear physics is dominant, it is difficult to make any 
strong argument for expecting frequency separation. In- 
deed, we would be surprised to not find violations of this 
assumption as we probe beyond the first few orders of 


4 



FIG. 3: Possible shapes from equations (3)-(4) for the mod- 
eled (3, 2) mode frequency [top panel] and amplitude [bottom 
panel] resulting from a merger with an at = 0.475 endpoint. 
All curves assume zero phase difference 5, and overall am- 
plitude is arbitrarily scaled to unity at t = 0. As with fig- 
ure 2, the two dashed (black) horizontal lines in the top panel 
mark the fundamental QNM frequencies wi (lower) and CJ32 
(higher) for a Kerr hole of the same final spin. 


magnitude in waveform precision. 

In the linear ringdown dynamics where this investiga- 
tion is focused, some degree of physical frequency sepa- 
ration can be expected, based on the separability of the 
Teukolsky equation, which describes small distortions of 
a stationary black-hole spacetime. The scale of physical 
linear mode-mixing can be quantified by careful consid- 
eration of quasinormal modes. 

The alternative hypothesis is that the mixing is an ar- 
tifact of our analysis, arising from choices that we make 
in setting up the angular mode decomposition. Perhaps 
our basis is not quite optimal, but we can find some other 
basis in which we more precisely recover frequency sep- 
aration. Indeed, given the freedom available in selecting 
such a representation, we have little grounds for suppos- 
ing that our first guess would be optimal. Here we con- 
sider two classes of choices in how to represent the space 
of gravitational radiation waveforms, which, in the full 


sense, has angular and retarded-time dimensions. 

The first choice we make is in how we define the 
spheres on which angular harmonic decomposition will 
be conducted. Within the structure of asymptotically 
flat spacetimes, gauge freedom in the choice of constant- 
retarded-time spheres can yield a frequency-dependent 
mode-mixing effect in the decomposed waveforms. This 
ambiguity arises from the freedom to re-parameterize the 
proper-time coordinate, the so-called “supertranslations” 
subgroup of the Bondi-Metzner-Sachs gauge group for 
outgoing radiation. We describe this possibility of super- 
translation gauge mixing in more detail below. We may 
generally expect that mode-mixing of this sort will be 
most evident in the late merger, where wavelengths are 
shortest. 

The next choice we make is in choosing the family of 
angular basis functions on the extraction spheres. In this 
case, the mixing arises if our chosen family of modal ba- 
sis functions used for radiation extraction differs from 
the optimal one in which frequency separation is best 
approximated. It is common to apply a spin-weighted 
spherical-harmonic basis, but a different choice may be 
motivated for the ringdown signals. Indeed, the separa- 
tion of the Teukolsky equation is not achieved in a spin- 
weighted spherical harmonic basis, but in a spin- weighted 
spheroidal- harmonic basis. Informally it is sometimes 
[BJK: cite for this?] assumed that this difference 
explains the sort of waveform phenomena we consider, 
though this has not been demonstrated. We label this 
effect angular-basis mixing. 

In the next subsections, we consider these possible mix- 
ing effects in detail, preparing for a quantitative study of 
the evidence for these effects in numerical data in Sec- 
tion V. 


A. Gauge Effects 

To understand the effect we are calling supertransla- 
tion gauge mixing, we must make a brief detour to de- 
scribe the gauge freedom in the representation of an out- 
going radiation field approaching future null infinity in an 
asymptotically flat spacetime. Consider such a spacetime 
in standard retarded-time coordinates { u , r, 9, cf>}. Scaled 
by r, the outgoing radiation field propagates outward on 
null rays labeled by u, 6 , and 4- Each polarization com- 
ponent can thus be described by a function of these vari- 
ables. The Bondi-Metzner-Sachs (BMS) [31, 32] group 
describes gauge transformations among these variables 
of the form 

6 ' = O' ( 6 , 4 ,), 4' = 4' (0,4), u' = K(0,4)(u- a(0,4)) , 

where (0,4>) — > (O' ,4') is a conformal transformation on 
a constant-M sphere with conformal factor K. 

For concreteness in the context of numerical relativity 
simulations, we note that it is common to make these 
gauge choices by specifying an “extraction sphere” lo- 
cated sufficiently far from the source where radiation field 
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calculations are realized. The effect of one class of BMS 
transformations, amounting to rotations of the extraction 
sphere, has been identified as an important concern when 
the choice of axis is not fixed by symmetry [33-36]. How- 
ever, the simulations in this study involve nonprecessing 
mergers, with no ambiguity in defining the orientation of 
the extraction sphere. 

But what happens if we make a small radial perturba- 
tion of the extraction sphere? It is clear that sufficiently 
small distortions of larger extraction spheres would have 
negligible impact on the intrinsic geometry of the sphere. 
The gauge effects of such distortions are described by a 
subset of the BMS transformations, known as supertrans- 
lations ■, with 9' = 9, <f>' = 0, and K = 1. 

Now consider the effect of a supertranslation on a grav- 
itational waveform tjj(u, 9, 0). Here we will make the ad- 
ditional assumption that a(9 , 0) is sufficiently small that 
we can approximate the effect of the supertranslation by 

0<y , 8, 0) ~ 0(«, 8, 0) + a(9, 0) ^-0(w, 9, 0), (5) 

and we can expand the supertranslation in terms of 
(scalar) spherical harmonics: 

a(0,0) = X>£M o^ M (M). (6) 

LM 

Then from (5), the measured radiation modes will be 
perturbed as follows: 

Q 

~ T 'y ^ (^)? (7) 

£.'m' 


where 


Cgrt 


= y ^pLM 


LM 


0 T l m -2 Y?' -lYr dtt 


= y ^pLM 


LM 


(2 L + 1) (2 1' + 1) 
47t (2/ + 1) 


1/2 


x (L, 0, l', 2| t, 2) (L, M, e, m'\e , m) . (8) 


In this paper we focus on mixing from the dominant 
mode, (£' = 2, m' = 2), with another m = 2 mode, fixing 
these values. For these cases the Clebsch-Gordan selec- 
tion rules require that M = 0 and l — 2 < L < t + 2. 
Then our mixing coefficient takes the form 


C ^222 — &lo 


’ 5 {2L + 1) 

_47t (21 + 1)_ 


1 1/2 


(L, 0,2,2|/,2) 2 . (9) 


For example, complete expansions for i = 3 and l = 4 
would yield 


C 3222 — Y — ^22%/3&io + 33-\/5&2 o + 22\/7&3o 

+22640 + \/Tff>50^ , 

C4222 = ^143"v/562o + 286"v/ 7630 + 702f>4o 

+9IVTT650 + 14vT3& 6 o) • 
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FIG. 4: Real and imaginary parts of p gaug e ,32 and p gauge ,42 
(for 620 = 1) for post-merger Kerr BHs of dimensionless spin 
at (and corresponding fundamental (2, 2) QNM frequency 
Mf(T22). 



The shape of the distorted extraction sphere is deter- 
mined by the coefficients &lo : for real a, we need the 
6 lo a lso to be real. The reality of the Clebsch-Gordan 
coefficients then implies that CV222 is also real. 

The other ingredient in the waveform-mode perturba- 
tion (7) is the derivative w.r.t. u on the right-hand side: 


= d u (A(«)e^“>) 





After merger, the effective coefficient yA/ A + iipj will 
asymptote to a constant complex number: 



1 

T’g'm' 


This implies a simple, QNM-driven leakage from the 
(2, 2) mode into higher-/ modes. Collecting terms, and 
working with the strain-rate h, during ringdown we have 

^(£,2) S ~ fyi,2) + iCi2 22 022 /l( 2,2) 

= h(e, 2 ) + Pgauge,« 2,2)- (10) 

In Fig. 4, we show the real and complex parts of the 
leakage parameters p gauge ,32 and p gau ge,42 for the sweep 
of end-state spins «f, assuming an unchanging scaling 
620 = 1 (and all other 6^ 0 = 0). The value of 620 is 
not physical, but gauge, and may differ between any two 
waveform determinations. The most important property 
we note is that the BMS leakage coefficients are nearly 
pure-imaginary at any fixed 620 and any spin Of. 
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B. Angular Basis Effects 

Another possible path to mixing arises from consider- 
ing what quasinormal-mode (QNM) frequencies actually 
represent. QNMs were discovered by Teukolsky in the 
context of the perturbation theory of Kerr black holes 
[37] . In developing this theory, Teukolsky worked with a 
background Kerr black hole in a very specific coordinate 
system due to Boyer & Lindquist [38]. 1 

A perturbed Kerr black hole will ring down to qui- 
escence through the emission of gravitational waves. 
These waves will have characteristic frequencies uu m 
and damping times T£ m given by the hole’s QNM spec- 
trum. 2 While the primary aim of QNM analysis is 
to determine the set of allowed complex frequencies 
(j£ m = ujdm + if Tim , these frequencies are tied to the 
radial and angular eigenfunctions arising from the sep- 
aration of the perturbation equations. These angular 
eigenfunctions are the spin-weighted spheroidal harmon- 
ics, - 2 y™{Mao; 9, <j)) = - 2 S™ (Mao; cos 9) e 1 ™^. 3 

Numerical waveform extraction from binary mergers, 
on the other hand, typically decomposes the waveforms 
onto the more generally motivated basis of spin- weighted 
spherical harmonics - 2 Y™ {9 , cf>) , which correspond to a 
spheroidal harmonic basis with Mao = 0: - 2 Y™($, cj>) = 
-2V™(0; 9, </>) [37]. Without an obvious nontrivial choice 
for Mao over the course of the evolving simulation which 
applies at all times, for all modes, decomposing with 
Mao — » 0 seems a natural choice. Here we consider an 
alternative choice, Mao — t Mf 0^022 • Using this basis re- 
quires knowing the final Kerr state (Mf , cq) of the merger 
before the decomposition can be applied, and the addi- 
tional task of numerically computing the basis functions 
(see Appendix A). Still this basis is not optimal for the 
subdominant modes. This unavoidable sub-optimality 
is discussed further in the next subsection. The distinc- 
tion between the spheroidal and spherical harmonics may 
be expected to yield the appearance of mode-mixing in 
the numerical waveform results even if we have elimi- 
nated the gauge freedom noted in the last section by op- 
timal correspondence with a suitably perturbed Boyer- 
Lindquist coordinate system. 

To estimate the apparent mode-mixing from this ba- 
sis mismatch, we can calculate the overlaps between the 
spheroidal harmonics (for a particular Mao) and the 
spherical harmonics. That is, we want to know the coef- 
ficients Si'im in 

OO 

2 yr (Mao- 9,(1 , ) = ^ 8 t>tm- 2 YF( 8 ,<l>). ( 11 ) 

i'=2 


1 Teukolsky theory can be reformulated on other backgrounds; see, 
e.g., [39]. 

2 We omit the principal quantum number n, assuming that we are 
dealing with the slowest-damped fundamental (n = 0) QNM. 

3 Here we use the symbol cr to denote a generic complex frequency. 
(7 Urn is a specific eigenvalue of the Kerr background. 


We describe our calculation of the - 2 $™ Ap- 
pendix A. To determine the overlaps we decompose 
the properly normalized spheroidal harmonic against the 
spherical harmonics in the usual way: 

svim =jdSl _ 2 yF(M { a i o 22 ; 9, </>) _ 2 0 )* 

= / dx - 2 SY l {Mia t o 2 2 ',x) - 2 S™(0;x)*. 

Now consider the idealized case where a physical ring- 
down signal is the simple combination of the fundamen- 
tal (2,2), (3,2), and (4,2) quasinormal modes (we omit 
Mao arguments for brevity): 

4 

h(t, r, 9, (j>) = ^ Hn{t,r ) _ 2 ^(M) 

£ 

4 

»X>'2(t,r)-2*?(M). (12) 

t' 

If we make the reasonable assumption that mixing i ^ £' 
products can be ignored for subdominant modes, then 
the measured spherical harmonic ringdown modes are ap- 
proximately: 



~ S 2 ' 22 % 22 (£, r), 


h^)(t,r) ; 

« Si'i 2 'Hi 2 (t, r) + Pba,sis,l2h 2 '2{t : 7"). 

(13) 

Here, the mixing coefficients are 



_ Si' 22 
Pbasis,^2 — 

(14) 


S2'22 

In Fig. 5, we plot the coefficients pbasis ,32 and Pbasis, 42 , 
evaluated at Mao = Mfaf 0 22 , where o 22 is the funda- 
mental QNM frequency of the (2, 2) mode for a Kerr hole 
of mass Mf and dimensionless spin Of. Note that (a) 
there is no ambiguity in overall scale for these coefficients 
(unlike the BMS-derived coefficients of the last section), 
and (b) they are strongly real-dominated. 

C. Physical mixing 

The discussion above exposes artifacts that arise from 
waveform decomposition using ordinary spin-weighted 
spherical harmonic functions. Here we ask whether, even 
with extraction spheres in the Boyer-Lindquist gauge, an- 
other decomposition using spin-weighted spheroidal har- 
monic functions can avoid mode-mixing. 

The question is non-trivial. Although each leading- 
order quasinormal ringdown mode exhibits angular de- 
pendence described by some kind of spin-weighted spher- 
ical harmonic angular function, they are not mutually 
given by the same kind of spin-weighted spherical har- 
monic angular functions, since each has its own distinct 
quasinormal frequency Oi m , and consequently a distinct 
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FIG. 5: Real and imaginary parts of the mixing coefficients 
/9basis,32 and pbasis.42 for post-merger Kerr BHs of dimension- 
less spin ac (and corresponding fundamental (2, 2) QNM fre- 
quency A7f (j 22 ) • At the zero-spin limit af — > 0, the leakage 
vanishes. 


preferred basis as labeled by Maa = Miotfcrtm- We must 
choose some particular orthonormal basis for the decom- 
position, and that basis cannot be simultaneously opti- 
mal for each mode. 

That the spheroidal harmonics associated with differ- 
ent QNM frequencies are not perfectly orthogonal has 
been demonstrated for high-spin Kerr holes by Berti et 
al. [40]. To quantify this for a general end-state spin af, 
we define new overlaps tin, between spheroidal harmon- 
ics associated with different m = 2 QNM frequencies: 

tin — dft -2T^(A7faf(j£2) -^yf (MfCkfcr^)*- (15) 

Figure 6 shows the magnitude of these overlaps for £ = 2 
(upper panel) and l = 3 (lower panel), and several val- 
ues of £'. From these plots, we see that the spheroidal 
harmonics for different Maa are not orthogonal, but 
show mixing by as much as ~ 4% for high spins (though 
the maximum overlaps occur at sub-maximal spins, as 
noted by [40]). The overlaps are also greatest for “near- 
est neighbor” modes: l = £' ± 1. For example, if we 
decomposed a waveform, including a non-trivial (2, 2) 
QNM, in the spheroidal basis corresponding to the (3, 2) 
mode ringdown frequency, then the corresponding curve 
in Fig. 6 would represent a mixing coefficient analogous 
to those in the previous subsections. There is no choice of 
orthonormal basis that will avoid all such mode mixing. 
In this sense, the angular non-orthogonality of the quasi- 
normal mode implies a form of physical mode-mixing, 
meaning that we can not perfectly isolate the QNM fre- 
quencies by any choice of angular basis. 

Fortunately it seems that the most evident mixing in- 
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FIG. 6: [Top panel] Magnitude of overlap te 2 between - 2 y£, 
evaluated at Maa = Mtata22 and - 2 TI, evaluated at 
Maa = Mfaian- [Bottom panel] Same for overlap te 3 be- 
tween -2T3 and -2 yl- By definition, |t32 1 = |i23 1 - 


volves the dominant (2, 2) mode frequency bleeding into 
higher-^ modes. With that assumption we may still elim- 
inate most physical mixing by choosing the basis compat- 
ible with this dominant quasinormal mode. If we decom- 
pose with the basis labeled by Maa = Mfafa 22 then 
the orthogonality of this particular basis will completely 
prevent the (2, 2) quasinormal mode from mixing in to 
any other decomposed modal waveform component. In 
this way we can eliminate any “physical mixing” of the 
particular form described in Section II. Mixing among 
subdominant modes, or mixing of subdominant modes 
into the decomposed (2, 2) waveform component will still 
occur at some level, but this is a smaller effect, which we 
do not focus on in this paper. 


IV. SIMULATIONS 

To investigate the mixing in a systematic way, we have 
surveyed several existing simulations of aligned-spin bi- 
naries, as well as carrying out new short simulations 


with the Goddard Hahndol evolution code. We choose 
our new black-hole binary (BHB) configurations in sev- 
eral groups of “merger-equivalent” classes, as described 
in B. The initial parameters for all these simulations, 
old and new, are presented in Table I. In Fig. 7, we 
show the distribution of these configurations as plots in 
the two-dimensional configuration-spaces {0:1,02} and 
{01,9}, where q = Mi/M 2 > 1 is the mass ratio, and 
a a = Sa/M\ is the dimensionless spin parameter of hole 
A, with physical values restricted to a a G [—1,1]. Many 
of the longer and higher-resolution evolutions have ap- 
peared in previous publications [25, 26]. Since our pri- 
mary interest here is strictly in the late-merger regime, 
newer evolutions begin only a few orbits before merger. 


A. Numerics 

The initial momenta of the newer evolutions were cho- 
sen by integrating the post-Newtonian equations of mo- 
tion, as outlined in [42, 43], with spin contributions to 
the Hamiltonian adapted from [44-48] , and the flux from 
[49]. Note that we did not attempt to reduce the eccen- 
tricity through tuning the initial momenta. 

The new evolutions use the Hahndol code paired with 
the “Curie” release of the Einstein Toolkit [50], incorpo- 
rating the Cactus Computational Toolkit [51] and the 
Carpet mesh-refinement driver [52]. 

In all cases, the initial data are of the standard Brandt- 
Briigmann puncture type [53] , using the Bowen- York [54] 
prescription for extrinsic curvature that exactly satis- 
fies the momentum constraint. We solve the remaining 
Hamiltonian constraint using the TwoPunctures spectral 
code [55]. 

To evolve these initial data, we employ the BSS- 
NOK 3+1 decomposition of Einstein’s vacuum equa- 
tions [56-58], with the alternative conformal variable 
W = e _2< ^ suggested in [59-61], constraint-damping 
terms suggested in [62], and the dissipation terms sug- 
gested in [63, 64], Our gauge conditions are the specific 
1+log lapse and Gamma-driver shift described in [65], 
which constitute a variant of the now-standard “moving 
punctures” approach [2, 3]. Our spatial derivatives use 
sixth-order-accurate differencing stencils, with the excep- 
tion of advection derivatives, where we use fifth-order- 
accurate mesh-adapted differencing (MAD) [66]. Our 
time- integration is performed with a fourth-order Runge- 
Kutta algorithm. 


B. Waveform Extraction 

We extract the gravitational waveforms from the sim- 
ulations through the radiative Weyl scalar ip 4 [28] . This 
is evaluated throughout the grid, and interpolated onto 
a set of coordinate spheres at extraction radii r G 
[40M, 90M] . Over each sphere, the interpolant is inte- 


grated against the set of spin- weighted spherical harmon- 
ics -2 Y™(9, (/>), up to £ = 5. 

In the extraction region, the grid spacing is between 
M/2 and 2 M, depending on the central resolution of the 
simulation. This is generally too coarse to resolve higher- 
frequency (and higher-m) modes with accuracy. Even for 
the dominant, relatively low-frequency, (2, ±2) modes, 
dissipation effects are visible that spoil the 1/r extrap- 
olation near and after merger. For this reason, we have 
used an r-extrapolation scheme that includes an explicit 
dissipative term in the amplitude of each mode: 

Aim{r) = o_ir + a 0 + a 2 r~ 2 , ^ m (r) = b 0 + b 2 r~ 2 . 

We have found this extrapolation procedure to be robust 
only for the higher-resolution simulations in this paper. 

As a result, a wavefonn-derived quantity / will have 
errors due to finite extraction radius and finite resolu- 
tion. For this paper, we make a very conservative error 
estimate by adding uncertainties linearly: 

A/ = A r f + A h /. 

For the finite-r error, we assume an uncertainty equal 
to the difference between the coefficient from the r- 
extrapolated highest-resolution data and that measured 
from the largest finite-r data at the same resolution. For 
finite-resolution error, we use the difference between the 
same-extraction-radius data at the coarse and fine reso- 
lutions as our estimate of the error in the fine-resolution 
result. For many configurations, we only have a single 
resolution available and the r-extrapolation is not reli- 
able at this resolution. For these, we adopt a conser- 
vative overall error estimate by taking the average error 
from comparable two- resolution configurations 4 and mul- 
tiplied it by 1.5. For amplitudes, this was a relative error, 
while for phase measurements, it was the absolute error. 


V. ANALYSIS OF WAVEFORMS 

Using the ringdown data from all the simulations in 
Section IV, we performed least-squares fits to the real 
part of the strain-rate (2,2) and (3,2) waveforms, using 
the forms of equations (l)-(2). Our fit is over the win- 
dow t G [20,55], where t = 0 is the time of peak (2,2) 
mode amplitude. By starting 20M after peak amplitude, 
we ensure that we are in the linear ringdown regime; 
by stopping at 55 M, we avoid the low-amplitude degra- 
dation seen in late-ringdown waveforms. As the tabu- 
lated version of the results would be excessively long, we 
present our raw results purely graphically. 


4 By “comparable” , we mean configurations that used the same 
numerical executable and grid structure, and whose lower- 
resolution version matched that of the single-resolution configu- 
rations. 
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TABLE I: Physical and numerical parameters of the initial data for all the runs presented, m i iP and roj, p are the bare puncture 
masses of the two pre-merger holes, ro is the initial coordinate separation, while Pot and Por are the initial transverse and 
radial components of the Bowen- York linear momentum. Madm is the total energy of the initial data, while the total infinite- 
separation mass of the system is estimated by the sum of the initial ADM masses of the individual holes [41] . We have found 
that for all cases here, this differs from the sum of apparent-horizon masses (calculated at times between t = 100 and 200), by 
less than a tenth of a percent. 


run name 

mi )P 

Wl2,p ^1 z 

S2z 

ro Pot 

(xlO 2 ) 

Por 

(xlO 4 ) 

Madm J2i M ADM,i 

Xl.UU 

0.301805 

0.301805 0.2 

0.2 

8.20 

10.32 

0.00 

0.988459 

1.000804 

Xl.mi 

0.454575 

0.454575 0.1 

0.1 

10.21 

9.25 

9.17 

0.99223 

1.002768 

XI .00 

0.487231 

0.487231 0.0 

0.0 

11.00 

9.01 

7.09 

0.990514 

1.000050 

X1JJD 

0.301805 

0.301805 0.2 

-0.2 

11.00 

9.01 

7.09 

0.990024 

0.999222 

XI. 5.00 

0.581359 

0.380645 0.0 

0.0 

7.12 

11.75 

29.17 

0.987252 

1.000000 

XI. 75.00 

0.619237 

0.345598 0.0 

0.0 

7.42 

11.01 

24.10 

0.988129 

1.000000 

X2.00 

0.649344 

0.314904 0.0 

0.0 

7.00 

11.00 

0.00 

0.987939 

1.000000 

X2JDU 

0.648662 

0.265507 -0.066666667 0.066666667 

10.00 

8.52 

7.63 

0.990951 

1.000009 

X2 . 5.00 

0.699349 

0.269501 0.0 

0.0 

7.40 

9.79 

20.53 

0.989664 

1.000000 

X3.00 

0.738687 

0.237505 0.0 

0.0 

8.88 

7.88 

8.96 

0.991673 

1.000000 

X4.00 

0.790000 

0.189000 0.0 

0.0 

8.47 

6.96 

0.00 

0.992912 

1.000310 

X5.U0 

0.822007 

0.157080 0.065083333 

0.0 

8.68 

5.91 

4.88 

0.993733 

1.000000 

X3_d0 

0.731667 

0.237705 -0.087566063 

0.0 

9.06 

7.84 

8.76 

0.99187 

1.000000 

X2_d0 

0.587677 

0.317821 -0.210380889 

0.0 

8.44 

9.93 

16.53 

0.989967 

1.000000 

X1.DD 

0.390411 

0.390411 -0.159125 

0.159125 

11.98 

8.84 

1.20 

0.990453 

0.998786 

X5.00 

0.824897 

0.157031 0.0 

0.0 

8.67 

5.97 

5.85 

0.993827 

1.000000 

X6.00 

0.848615 

0.133064 0.0 

0.0 

7.55 

5.84 

6.94219 

0.994008 

1.000000 

X5JD0 

0.822405 

0.156318 -0.052232639 

0.0 

8.09 

6.32 

7.00085 

0.993556 

1.000000 

X4.D0 

0.778549 

0.188766 -0.1213184 

0.0 

8.57 

7.04 

7.78076 

0.992926 

1.000000 

X3JD0 

0.692530 

0.237756 -0.21614625 

0.0 

9.17 

7.93 

8.80172 

0.99219 

1.000000 

X2JDD 

0.531347 

0.260245 -0.277766667 -0.069441667 

10.72 

8.56 

7.81844 

0.992008 

1.000000 
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FIG. 7: Simulated configurations from Table I, represented as points in two-dimensional {ai, « 2 }-space [left panel] and (ai, q}- 
space [right panel]. Note that some points are associated with multiple simulations. 
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We begin by showing the nature of the complex numer- 
ical “leakage parameter” derived from the ratio of fitted 
parameters from the measured (2, 2) and (3, 2) modes 
during ringdown, using (l)-(2): 


= p 3 2A 22 e i ^ t+s ^ 

A 22 e^ a22t+S2 A 


(16) 


Figure 8 shows the real and imaginary parts of this leak- 
age for all configurations presented in this paper, as a 
function of the dimensionless spin af of the post-merger 
hole. 


A. Comparing Hypotheses 

In Section 111 we discussed two possible causes for 
mode-mixing effects of the form 

= A e2 e l( ' ae2t+Se2) + pe2A 2 2e i(a22t+S22) , (17) 

described in Section II. If the mixing is caused by BMS 
supertranslation gauge ambiguity, then we would expect 
nearly pure imaginary pi 2 . On the other hand, if the 
mixing derives from the distinction between spheroidal 
and spherical harmonic angular functions, then we expect 
predominantly real pi 2 of a quantified size. In Fig. 8 we 
see that the argument of p nU m,32 is close to zero, within 
error bars for most cases, making p n um,32 predominantly 
real, consistent with the spheroidal harmonic hypothesis. 
The largest deviations from zero are also those with the 
largest uncertainties arising from the QNM fit process. 

The analysis of Section III A suggests that a change of 
supertranslation gauge would give rise to mode-mixing 
coefficients with a numerically significant imaginary part 
in the measured waveform. Since the imaginary part of 
Pnum,32 is so small, any supertranslation gauge effects are 
negligible at the level of interest here. We can estimate 
the degree of gauge constraint implied by a null mea- 
surement of this effect. From the bottom panel of Fig. 8, 
one sees that the imaginary component of the mixing co- 
efficient Pnum ,32 is constrained to values within ±0.02 in 
almost all cases. If we generously assumed that all of this 
imaginary mixing was caused by gauge distortion of the 
extraction sphere, by comparison with Fig. 4, we would 
conclude that the amplitude of the distortion (b 2 o specif- 
ically) would be have to be smaller than about 0.2 M, 
suggesting a remarkable level of supertranslation gauge 
optimality in these simulations. 



FIG. 8: Real [top panel] and imaginary [bottom panel] parts 
of complex leakage parameter p n um ,32 (16), shown as a func- 
tion of ctf, the dimensionless spin of the post-merger hole for 
all runs, indicating that in almost all cases, the leakage pa- 
rameter is predominantly real. 


with the magnitude of the leakage coefficients Pbasis ,32 
(14) plotted in Fig. 5 (blue solid curve). From the close 
fit, it appears that the leakage is in fact dominated by 
this spheroidal/spherical harmonic mismatch. That is, 
even though the post-merger background coordinate sys- 
tem should not be expected to closely resemble the Kerr- 
Boyer-Lindquist slicing assumed by Teukolsky’s pertur- 
bative work, nevertheless, this expected warping is not 
as important as our choice of harmonic basis functions. 


B. Testing the Spheroidal Leakage Model 

We have seen that the numerical results for the com- 
plex argument of p% 2 are consistent with the spherical- 
spheroidal mixing hypothesis, but this hypothesis also 
makes quantitative predictions for the magnitude I/O 32 1 - 
In the top panel of Fig. 9, we plot the magnitude of 
Pnum,32 as a function of af. We overlay these points 


The bottom panel of Fig. 9 shows the complex ampli- 
tude of the equivalent parameter p numj 42 governing the 
leakage of the (2, 2) mode into the measured (4, 2) mode. 
Although this is also consistent with expectations from 
angular-basis mixing (blue solid curve), the relative er- 
rors swamp the numerical data, and higher-resolution nu- 
merics will be needed to establish the relation unambigu- 
ously. 
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FIG. 9: [Top panel] Amplitude of measured (2, 2) — t (3, 2) 
leakage parameter p n um ,32 as a function of at, the dimension- 
less spin of the post-merger hole for all runs. The curve is the 
theoretical ratio Pbasis.32 due to spheroidal-spherical leakage. 
[Bottom panel] Same for (2, 2) — » (4, 2) leakage ratio p n um, 42 - 


C. Finding the Residual (3,2) Mode Amplitude 

If we regard the measured (3,2) mode as the combi- 
nation of a “true” (3,2) mode A32 exp 7(032! + fe) and 
a piece of the (2, 2) mode, we may ask whether we can 
model the residual (3,2) contribution. When looking at 
the entire suite of simulations, it is difficult to see a dis- 
tinct pattern in these true (3, 2) amplitudes. However, it 
is instructive to carry out a particular slice in configura- 
tion space. 

In Fig. 10, we show a subset of the (3, 2) amplitudes 
formed by the mergers of nonspinning binaries, with 
mass ratio q = Mi/M 2 £ {1.0, 6.0}. Error bars in 
this plot have been estimated in the same way as for 
Fig. 9. Clearly the high-g behavior seems to decay to 
some constant amplitude, while there is some local min- 
imum around 77 = 0.21 (between q = 2 and q = 2.25), in- 
dicating that perhaps at this mass-ratio, the (3, 2) QNM 
is hardly excited at all. 



FIG. 10: Corrected amplitude A 32 of the (3, 2) strain-rate 
mode for nonspinning binaries with mass ratio q £ {1.0, 6.0}. 
The solid (green) curve is a fit to these points of the functional 
form (18). 

We also present an empirical fit to this data of the 
functional form 

^32(77) = yj (a — be~ x / r ') 2 + c 2 , (18) 

where the parameters take the values a = 0.0147±0.0002, 
b = 1.5 ± 0.1, c = 0.0026 ± 0.0002, and A = 0.98 ± 0.02. 
Fits of this form are expected to be useful for generating 
merger template waveforms for the subdominant modes. 

VI. DISCUSSION 

In this paper, we have investigated “bumps” measured 
in the merger-ringdown portion of certain gravitational- 
radiation angular waveform modes from the numeri- 
cal simulation of the coalescence of black-hole binaries 
(BHBs). These bumpy modes appear to contain signifi- 
cant contributions from the dominant (2, 2) mode, indi- 
cating some kind of mode-mixing at work. 

We have considered three classes of effects that may 
contribute to mode-mixing in numerically extracted and 
decomposed merger-ringdown waveforms. These are: 
gauge effects, arising from supertranslation gauge free- 
dom for outgoing radiation in general asymptotically 
flat spacetimes (see Sec. Ill A); angular-basis effects, re- 
lating to a choice between spin-weighted spherical or 
quasi-normal-mode-adapted spheroidal harmonic bases 
(Sec. IIIB); and physical quasinormal-mode mixing ef- 
fects that are independent of any representation changes 
(Sec. IIIC). 

We have identified and analyzed the measured mode- 
mixing bumps in the most prominent of the bumpy grav- 
itational waveform modes modes — £ = 3,m = 2 — mea- 
sured from a set of numerical evolutions of aligned-spin 
BHB mergers. Our analysis has allowed us to distinguish 
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between the contributions of our three mode-mixing ef- 
fects. We find that the angular-basis effects dominate. 
Although other kinds of effects may be present - like the 
frequency-dependent gauge supertranslations discussed 
in Section III A - they cannot be seen clearly here with 
the level of accuracy available from our current simula- 
tions. 

In this way our analysis further codifies the results 
from the ringdown stage of the aligned-spin mergers. 
This was originally prompted by our work on a multi- 
mode waveform model based on the implicit rotating 
source (IRS) picture of black- hole merger [25, 26]. In this 
model, the dominant and leading sub-dominant wave- 
form modes from binary mergers were seen to share a 
common rotational phase, with a corresponding rota- 
tional frequency that increased monotonically through 
inspiral and merger, plateauing during ringdown. The 
corresponding mode amplitudes could be modeled by a 
simple, few-parameter functional form that depends on 
the frequency function, with a single well-defined peak. 
Attempting to extend this to the (3, 2) mode proved prob- 
lematic, as the measured mode was no longer monotonic 
in frequency, or single-peaked in amplitude. 

More broadly, we expect our results to provide guid- 
ance in the ongoing effort of combining results of ana- 
lytic and numerical relativity studies toward the goal of a 
fully developed family of efficient and accurate black-hole 
merger waveforms. Because the comparison of waveform 
models is typically conducted mode by mode in decom- 
posed form, the issues we have studied may lead to un- 
necessarily spurious features in particular waveform rep- 
resentations. 

We estimate, for instance, that supertranslation gauge 
changes that would effectively distort the shape of ar- 
bitrarily large waveform-extraction spheres on scales of 
order Mf or smaller would be sufficient to qualitatively 
influence the mode-mixing features focused on in this 
study. The absence of such effects is itself intriguing, 
suggesting that we have achieved nearly optimal choice 
of supertranslation gauge. Our near-optimal spheroidal 
harmonic basis is consistent with quasinormal-mode dis- 
tortions of Kerr space-time in the Boyer-Lindquist co- 
ordinate system. That we see negligible supertransla- 
tion mode-mixing suggests that the outer regions of our 
numerical space-times asymptotically approach distorted 
Kerr in Boyer-Lindquist coordinates faster (in powers 
of 1/r) than the asymptotic approach to perturbed 
Minkowski spacetime. This seems plausible, based on 
our choice of numerical gauge, which approximates max- 
imal time slicing and T* = 0 spatial coordinates. The 
latter condition will yield spatially isotropic coordinates 
where possible. 

Nonetheless, it seems that we have been lucky to stum- 
ble onto a near-optimal representation as other incompat- 
ible gauge choices may also be reasonable in the numer- 
ical simulation context. In continued pursuit of higher- 
precision waveform comparisons and higher-fidelity ana- 
lytic models (see, e.g., the NR-AR project [67]), we ex- 


pect such considerations to grow in significance. (They 
may also be crucial in studies of how the pre-merger 
BHB configuration is encoded in the relative amplitude 
of different quasi-normal modes during ringdown; see, 
e.g. [68].) Similarly we find that physical mode-mixing 
among the quasinormal modes will prevent any ortlionor- 
mal representation from fully separating frequencies at 
sufficiently high precision. 
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Appendix A: Calculating Spheroidal Harmonics 

As there are no closed- form solutions for the -23^™; 
we must proceed numerically. While setting up the pop- 
ular continued-fraction method for computing the QNM 
frequencies of a Kerr BH, Leaver [69] presents the follow- 
ing power-series expansion for the polar-angle function 
_2.S7 1 (M aa; cos 9), due originally to Baber & Hasse [70] 
(we specialize here to s = — 2): 

2 SF(Maa ; a:) =e Maax {l + a;) |m+2|/2 (l - ;c )l— 2 I/ 2 

OO 

x^a n (l + x) n , (Al) 

n = 0 

where the expansion coefficients a n are determined up 
to an overall scaling — the value of ao — by the same 
recurrence relations that yield the QNM frequencies. For 
our desired Kerr spin at, we first determine the (complex) 
fundamental QNM frequency of the (2,2) mode, Mf<722- 
Next, assuming ao = 1, we use the recurrence relations 
from [69] to determine the a n (in practice, we truncate 
the series at n = 14). Requiring that 

J dx\ - 2 Sy l (MfafCT 22 ', x)\ 2 = 1 

then fixes ao, supplying the correct normalization of the 

«n • 

Appendix B: Kerr-Equivalent Black-Hole Binaries 

The end-point of any merger of BHBs in vacuum is 
expected to be a single Kerr black hole, parameterized by 
two numbers, the mass Mf and spin angular momentum 
Sf = afM f 2 . These should satisfy the global conservation 
rules: 

Mf = MadM — -Erad, (Bl) 

= JaDM 3r ac f, (B2) 


13 


where Madm and Jadm are the Arnowitt-Deser-Misner 
(ADM) energy and total angular momentum of the ini- 
tial data, and E la d and J ra d are the energy and angular 
momentum emitted in gravitational radiation during the 
course of the evolution. 

Fixing the initial separation of the binary, and taking 
its total mass to be M = Mi + M 2 (> Madm for any 
finite initial separation), and assuming zero eccentricity, 
the black-hole binary will have seven free parameters: 
{g, Si, £2}, where q = Mi/M 2 > 1 is the mass ratio, 
and Sa are the spin angular momentum vectors of the 
two holes. However, the end-state has just two param- 
eters: {Mf,5f}, so there must be a large degeneracy in 
the initial parameters. 

Viewing the BHB coalescence as a kind of simple par- 
ticle interaction, Boyle et al. [71] used symmetry argu- 
ments to restrict the possible end-states of the BHB 
merger. This is the basis of end-state formulae by Tichy 
& Marronetti [72]. Other models have been developed 
by Buonanno et al. [73], Lousto et al. [74], Barausse & 
Rezzolla [75, 76] and others. 

In the case of initially orbit-aligned spins, the initial 
parameter space is three-dimensional: {g, S\,S 2 }. We 
use the simplest applicable formulas for the achieved end- 
state for an aligned-spin system. The end-state mass 
formula we take from Eq. (5) of [74]: 

Mf = 1 — qEisco ~ E 2 ?f — E 3 rf 
V 1 2 

~ 2 [ E s( a 2 + 9 2 4 «i) + E S {1 ~ q){a 2 - gap) 

+E A {a 2 + qai ) 2 + E D {a 2 - gai) 2 ] , (B3) 

where 17 = MjM 2 /(Mi+M 2 ) 2 = g/(l+g) 2 is the symmet- 
ric mass ratio of the binary, and the fitting parameters 
are: 


Fusco =1 


y/8 

3 


+ 0.103803?7 + 


(q(l + 2g)aq + (2 + g)a 2 ) 

36^3(1 + g) 2 


5(g«i - q 2 ) 2 
+ 162^/2(1 + g) 2 ’ 

E 2 =0.341, E 3 = 0.522, E s = 0.673, 

E s = - 0.36, E a = -0.014, E D = 0.26. 


For the final spin, one model with just enough com- 


plexity for our data sets here was given by [75, 76] 5 : 

af = a+s^qa 2 +S5T] 2 a+toria+2\/3ri+t 2 ri 2 + t 3 r/ 3 , (B4) 

where a = (g 2 aq + a 2 )/(q 2 + 1) and the coefficients 
{+4? 85, fp, f 2 , £3} are. 


s 4 = -0.1229+ 0.0075, s 5 = 0.4537+ 0.1463, 
t 0 = -2.8904 ± 0.0359, t 2 = -3.5171 ± 0.1210, 
t 3 = 2.5763 + 0.4833. 

TABLE II: Final mass and spin of the post-merger Kerr BH, 
as measured by radiation balance (Mf, at), and as predicted 
by phenomenological equations (B3)-(B4) (MqRiTjaf.AEi)- 
The final two columns give the percentage relative error be- 
tween the measured and predicted values, which never exceeds 
1.6 % for the mass and 2.1 % for the spin. 


run name Mf oif 

AffifjjT <*f,AEI 

SM t (%) Sat (%) 

Xl.UU 0.9156 0.9053 

0.9287 0.9112 

1.43 0.65 

Xl_uu 0.9393 0.8119 

0.9391 0.8038 

0.03 0.99 

X1.00 0.9520 0.6886 
Xl.UD 0.9505 0.6839 

0.9497 0.6865 
0.9359 0.6865 

0.24 0.31 

1.54 0.38 

XI. 5.00 0.9558 0.6664 

0.9534 0.6644 

0.25 0.30 

XI. 75.00 0.9588 0.6475 

0.9565 0.6452 

0.24 0.35 

X2.00 0.9614 0.6254 
X2.DU 0.9610 0.6120 

0.9596 0.6244 
0.9559 0.6244 

0.19 0.17 

0.54 2.02 

X2.5.00 0.9671 0.5833 

0.9654 0.5824 

0.18 0.16 

X3.00 0.9716 0.5432 

0.9702 0.5429 

0.15 0.07 

X4.00 0.9782 0.4780 
X5.U0 0.9816 0.4741 
X3J)0 0.9737 0.4735 
X2.D0 0.9683 0.4704 
Xl.dd 0.9646 0.4825 

0.9812 0.4748 
0.9773 0.4748 
0.9720 0.4760 
0.9649 0.4765 
0.9674 0.4786 

0.31 0.68 

0.44 0.15 

0.18 0.52 

0.35 1.31 

0.30 0.81 

X5.00 0.9826 0.4186 

0.9821 0.4202 

0.06 0.37 

X6.00 0.9857 0.3718 
X5.D0 0.9834 0.3736 
X4_D0 0.9803 0.3728 
X3_D0 0.9762 0.3697 
X2_DD 0.9718 0.3788 

0.9854 0.3762 
0.9791 0.3762 
0.9791 0.3762 
0.9739 0.3762 
0.9729 0.3762 

0.02 1.18 

0.44 0.68 

0.13 0.91 

0.23 1.75 

0.11 0.71 


Using these formulae, we have constructed a set of con- 
figurations, which we present in Table I, grouped by final 
Kerr spin. 


[1] F. Pretorius, Phys. Rev. Lett. 95, 121101 (2005), 
arXiv:gr-qc/0507014. 

[2] M. Campanelli, C. O. Lousto, P. Marronetti, and Y. Zlo- 
chower, Phys. Rev. Lett. 96 , 111101 (2006), arXiv:gr- 
qc/0511048. 

[3] J. G. Baker, J. M. Centrella, D.-I. Choi, M. Koppitz, 
and J. R. van Meter, Phys. Rev. Lett. 96 , 111102 (2006), 
arXiv:gr-qc/0511103. 

[4] J. G. Baker, M. Campanelli, F. Pretorius, and Y. Zlo- 


chower, Class. Quantum Grav. 24 . S25 (2007), arXiv:gr- 
qc/0701016. 

[5] M. D. Hannam, S. Husa, J. G. Baker, M. Boyle, 
B. Briigmann, T. Chu, E. N. Dorband, F. Herrmann, 
I. Hinder, B. J. Kelly, et al., Phys. Rev. D 79 . 084025 
(2009), arXiv:0901.2437 [gr-qc]. 

[6] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 
Rev. D 74, 041501(R) (2006), arXiv:gr-qc/0604012. 

[7] A. Buonanno, G. B. Cook, and F. Pretorius, Phys. Rev. 


14 


D 75, 124018 (2007), arXiv:gr-qc/0610122. 

[8] M. Boyle, D. A. Brown, L. E. Kidder, A. H. Mroue, H. P. 
Pfeiffer, M. A. Scheel, G. B. Cook, and S. A. Teukolsky, 
Phys. Rev. D 76, 124038 (2007), arXiv:0710.0158 [gr-qc], 

[9] M. A. Scheel, M. Boyle, T. Chu, L. E. Kidder, K. D. 
Matthews, and H. P. Pfeiffer, Phys. Rev. D 79, 024003 
(2009), arXiv:0810.1767 [gr-qc]. 

[10] J. A. Gonzalez, U. Sperhake, and B. Briigmann, Phys. 
Rev. D 79, 124006 (2009), arXiv:0811.3952 [gr-qc]. 

[11] C. O. Lousto, H. Nakano, Y. Zlochower, and 
M. Campanelli, Phys. Rev. Lett. 104, 211101 (2010), 
arXiv:1001.2316 [gr-qc]. 

[12] J. M. Centrella, J. G. Baker, B. J. Kelly, and 
J. R. van Meter, Rev. Mod. Phys. 82, 3069 (2010), 
arXiv:1010.5260 [gr-qc]. 

[13] I. Hinder, Class. Quantum Grav. 27, 114004 (2010), in- 
vited paper from Numerical Relativity and Data Anal- 
ysis (NRDA) 2009, Albert Einstein Institute, Potsdam, 
arXiv:1001.5161 [gr-qc]. 

[14] P. Ajith, M. D. Hannam, S. Husa, Y. Chen, 

B. Briigmann, E. N. Dorband, D. Muller, F. Ohme, 

D. Pollney, C. Reisswig, et al., Phys. Rev. Lett. 106, 

241101 (2011), arXiv:0909.2867 [gr-qc]. 

[15] Y. Pan, A. Buonanno, L. T. Buchman, T. Chu, L. E. 
Kidder, H. P. Pfeiffer, and M. A. Scheel, Phys. Rev. D 
81, 084041 (2010), arXiv:0912.3466 [gr-qc]. 

[16] L. Santamaria, F. Ohme, P. Ajith, B. Briigmann, E. N. 
Dorband, M. D. Hannam, S. Husa, P. Mosta, D. Poll- 
ney, C. Reisswig, et al., Phys. Rev. D 82, 064016 (2010), 
arXiv: 1005. 3306 [gr-qc]. 

[17] K. G. Arun, B. R. Iyer, B. S. Sathyaprakash, 

and S. Sinha, Phys. Rev. D 75, 124002 (2007), 
arXiv:0704.1086 [gr-qc]. 

[18] K. G. Arun, B. R. Iyer, B. S. Sathyaprakash, S. Sinha, 
and C. Van Den Broeck, Phys. Rev. D 76, 104016 (2007), 
arXiv:0707.3920 [astro-ph], 

[19] M. Trias and A. M. Sintes, Class. Quantum Grav. 25, 
184032 (2008), 12th GWDAW (Gravitational Wave Data 
Analysis Workshop), arXiv:0804.0492 [gr-qc]. 

[20] S. T. McWilliams, J. I. Thorpe, J. G. Baker, and B. J. 
Kelly, Phys. Rev. D 81, 064014 (2010), arXiv:0911.1078 
[gr-qc]. 

[21] J. N. Goldberg, A. J. MacFarlane, E. T. Newman, 
F. Rohrlich, and E. C. G. Sudarshan, J. Math. Phys. 
8, 2155 (1967). 

[22] Y. Wiaux, L. Jacques, and P. Vandergheynst, J. Comp. 
Phys. 226, 2359 (2007), arXiv:astro-ph/0508514. 

[23] E. Berti, V. Cardoso, J. A. Gonzalez, U. Sperhake, M. D. 
Hannam, S. Husa, and B. Briigmann, Phys. Rev. D 76, 
064034 (2007), arXiv:gr-qc/0703053. 

[24] E. Berti, V. Cardoso, J. A. Gonzalez, U. Sperhake, and 
B. Briigmann, Class. Quantum Grav. 25, 114035 (2008), 
arXiv:071 1.1097 [gr-qc]. 

[25] J. G. Baker, W. D. Boggs, J. M. Centrella, B. J. Kelly, 
S. T. McWilliams, and J. R. van Meter, Phys. Rev. D 
78, 044046 (2008), arXiv:0805.1428 [gr-qc]. 

[26] B. J. Kelly, J. G. Baker, W. D. Boggs, S. T. McWilliams, 
and J. M. Centrella, Phys. Rev. D 84, 084009 (2011), 
arXiv: 1107. 1181 [gr-qc]. 

[27] Y. Pan, A. Buonanno, M. Boyle, L. T. Buchman, L. E. 
Kidder, H. P. Pfeiffer, and M. A. Scheel, Phys. Rev. D 
84, 124052 (2011), arXiv: 1106.1021 [gr-qc]. 

[28] J. G. Baker, M. Campanelli, and C. O. Lousto, Phys. 
Rev. D 65, 044001 (2002), arXiv:gr-qc/0104063. 


[29] J. G. Baker, M. Campanelli, C. O. Lousto, and R. Taka- 
hashi, Phys. Rev. D 65, 124012 (2002), arXiv:astro- 
ph/0202469. 

[30] J. G. Baker, J. M. Centrella, D.-I. Choi, M. Koppitz, 
and J. R. van Meter, Phys. Rev. D 73, 104002 (2006), 
arXiv:gr-qc / 0602026. 

[31] H. Bondi, M. G. J. van der Burg, and A. W. K. Metzner, 
Proc. R. Soc. London Ser. A 269, 21 (1962). 

[32] R. K. Sachs, Proc. R. Soc. London Ser. A 270, 103 
(1962). 

[33] P. Schmidt, M. D. Hannam, S. Husa, and P. Ajith, Phys. 
Rev. D 84, 024046 (2011), arXiv:1012.2879 [gr-qc], 

[34] R. O’Shaughnessy, B. Vaishnav, J. Healy, Z. Meeks, and 
D. M. Shoemaker, Phys. Rev. D 84, 124002 (2011), 
arXiv:1109.5224 [gr-qc]. 

[35] R. O’Shaughnessy, J. Healy, L. London, Z. Meeks, and 
D. M. Shoemaker, Phys. Rev. D 85, 084003 (2012), 
arXiv:1201.2113 [gr-qc]. 

[36] P. Schmidt, M. D. Hannam, and S. Husa (2012), 
arXiv: 1207.3088 [gr-qc]. 

[37] S. A. Teukolsky, Phys. Rev. Lett. 29, 1114 (1972). 

[38] R. H. Boyer and R. W. Lindquist, J. Math. Phys. 8, 265 
(1967). 

[39] M. Campanelli, G. Khanna, P. Laguna, J. A. Pullin, 
and M. P. Ryan, Class. Quantum Grav. 18, 1543 (2001), 
arXiv:gr-qc/0010034. 

[40] E. Berti, V. Cardoso, and M. Casals, Phys. Rev. D 73, 
024013 (2006), Erratum: ibid. 73, 109902(E) (2006), 
arXiv:gr-qc/0511111. 

[41] B. D. Baker (2002), arXiv:gr-qc/0205082. 

[42] S. Husa, M. D. Hannam, J. A. Gonzalez, U. Sperhake, 
and B. Briigmann, Phys. Rev. D 77, 044037 (2008), 
arXiv:0706.0904 [gr-qc]. 

[43] M. Campanelli, C. O. Lousto, H. Nakano, and Y. Zlo- 
chower, Phys. Rev. D 79, 084010 (2009), arXiv:0808.0713 
[gr-qc]. 

[44] A. Buonanno, Y. Chen, and T. Damour, Phys. Rev. D 
74, 104005 (2006), arXiv:gr-qc/0508067. 

[45] T. Damour, P. Jaranowski, and G. Schafer, Phys. Rev. 
D 77, 064032 (2008), arXiv:0711.1048 [gr-qc]. 

[46] R. A. Porto and I. Z. R.othstein, Phys. Rev. Lett. 97, 
021101 (2006), arXiv:gr-qc/0604099. 

[47] J. Steinhoff, S. Hergt, and G. Schafer, Phys. Rev. D 77, 
081501(R) (2008), arXiv:0712.1716 [gr-qc], 

[48] J. Steinhoff, S. Hergt, and G. Schafer, Phys. Rev. D 78, 
101503(R) (2008), arXiv:0809.2200 [gr-qc], 

[49] L. Blanchet, A. Buonanno, and G. Faye, Phys. Rev. D 
74, 104034 (2006), Errata: ibid. 75, 049903(E) (2007); 
81, 089901(E) (2010), arXiv:gr-qc/0605140. 

[50] F. Lofffer, J. Faber, E. Bentivegna, T. Bode, P. Diener, 
R. Haas, I. Hinder, B. C. Mundim, C. D. Ott, E. Schnet- 
ter, et al., Class. Quantum Grav. 29, 115001 (2011), 
arXiv: 11 11.3344 [gr-qc]. 

[51] Cactus computational toolkit home page , http : //www . 
cactuscode . org. 

[52] Carpet: Adaptive Mesh Refinement for the Cactus 

Framework, http : / /www . carpetcode . org/. 

[53] S. R. Brandt and B. Briigmann, Phys. Rev. Lett. 78, 
3606 (1997), arXiv:gr-qc/9703066. 

[54] J. M. Bowen and J. W. York Jr., Phys. Rev. D 21, 2047 
(1980). 

[55] M. Ansorg, B. Briigmann, and W. Tichy, Phys. Rev. D 
70, 064011 (2004), arXiv:gr-qc/0404056. 

[56] T. Nakamura, K.-I. Oohara, and Y. Kojima, Prog. Theor. 


15 


Phys. Suppl. 90, 1 (1987). 

[57] M. Shibata and T. Nakamura, Phys. Rev. D 52, 5428 
(1995). 

[58] T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D 59, 
024007 (1998), arXiv:gr-qc/9810065. 

[59] J. R. van Meter, in From Geometry to Numerics (Institut 
Henri Poincare, Paris, 2006), http://luth2.obspm.fr/ 
IHP06/workshops/ geomnum/ si ides/ vanmeter . pdf . 

[60] W. Tichy and P. Marronetti, Phys. Rev. D 76, 061502(R.) 
(2007), arXiv:gr-qc/0703075. 

[61] P. Marronetti, W. Tichy, B. Briigmann, J. A. Gonzalez, 
and U. Sperhake, Phys. Rev. D 77, 064010 (2008), 
arXiv : 0 709 . 2 1 60 [gr-qc] . 

[62] M. D. Duez, S. L. Shapiro, and H.-J. Yo, Phys. Rev. D 
69, 104016 (2004), arXiv:gr-qc/0401076. 

[63] H.-O. Kreiss and J. Oliger, Methods for the approximate 
solution of time dependent problems, no. 10 in GARP 
Publications Series (World Meteorological Organization 
and International Council of Scientific Unions, Geneva, 
1973). 

[64] P. Hiibner, Class. Quantum Grav. 16, 2823 (1999), 
arXiv:gr-qc/9903088. 

[65] J. R. van Meter, J. G. Baker, M. Koppitz, and D.-I. Choi, 
Phys. Rev. D 73, 124011 (2006), arXiv:gr-qc/0605030. 

[66] J. G. Baker and J. R. van Meter, Phys. Rev. D 72, 104010 
(2005), arXiv:gr-qc/0505100. 


[67] The numerical relativity and analytical relativity (NRAR) 
collaboration, https : / / www . ninja-project . org/ doku. 
php? id=nr ar : home . 

[68] I. Kamaretsos, M. D. Hannam, and B. S. Sathyaprakash, 
Phys. Rev. Lett. 109, 141102 (2012), arXiv:1207.0399 
[gr-qc]. 

[69] E. W. Leaver, Proc. R. Soc. London Ser. A 402, 285 
(1986). 

[70] W. G. Baber and H. R. Hasse, Mathematical Proceedings 
of the Cambridge Philosophical Society 31, 564 (1935). 

[71] L. Boyle and M. Kesden, Phys. Rev. D 78, 024017 (2008), 
arXiv:0712.2819 [astro-ph], 

[72] W. Tichy and P. Marronetti, Phys. Rev. D 78, 081501(R) 
(2008), arXiv:0807.2985 [gr-qc], 

[73] A. Buonanno, L. E. Kidder, and L. Lehner, Phys. Rev. 
D 77, 026004 (2008), arXiv:0709.3839 [astro-ph], 

[74] C. O. Lousto, M. Campanelli, Y. Zlochower, and 
H. Nakano, Class. Quantum Grav. 27, 114006 (2010), 
arXiv:0904.3541 [gr-qc]. 

[75] L. Rezzolla, E. Barausse, E. N. Dorband, D. Pollney, 
C. Reisswig, J. Seiler, and S. Husa, Phys. Rev. D 78, 
044002 (2008), arXiv:0712.3541 [gr-qc]. 

[76] E. Barausse and L. Rezzolla, Astrophys. J. 704, L40 
(2009), arXiv:0904.2577 [gr-qc]. 


